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Abstract 

We present the results of the calculation of the 'photon' + 1 jet rate at (D(aa s ) for 
LEP energies. By comparing these results with the data from the ALEPH Collabora- 
tion we make a next-to-leading order determination of the quark-to-photon fragmenta- 
tion function D g ^ 7 (z,/ip) at 0{aa s ). The predictions obtained using this fragmenta- 
tion function for the isolated rate, defined as the 'photon' + 1 jet rate for z > 0.95, are 
found in good agreement with the ALEPH data. The next-to-leading order corrections 
are moderate demonstrating the perturbative stability of this particular isolated pho- 
ton definition. We have also computed the inclusive photon energy distribution and 
found good agreement with the OPAL data. 



*Talk presented at the Zeuthen Workshop on Elementary Particle Theory "Loops and Legs in Gauge 
Theories" Rheinsberg, Germany, April 19-24, 1998. 



1 Introduction 



Photons produced in hadronic collisions arise essentially from two different sources: the di- 
rect production of a photon off a primary parton or through the fragmentation of a hadronic 
jet into a single photon carrying a large fraction of the jet energy. The former gives rise to 
perturbatively calculable short-distance contributions whereas the latter is primarily a long 
distance process which cannot be calculated within perturbative QCD. It is described by 
the process-independent parton-to-photon fragmentation function |L| which must be deter- 
mined from experimental data. Its evolution with the factorization scale fip can however be 
determined by perturbative methods. 

The ALEPH Collaboration at CERN has measured the quark-to-photon fragmentation 
function D q ^ || from an analysis of two jet events in which one of the jets contains a 
photon carrying a large fraction (z > 0.7) of the jet energy. These 'photon' +1 jet events are 
defined by the application of the Durham jet clustering algorithm |j| to both the hadronic 
and electromagnetic clusters. In this democratic approach, the photon is called isolated 
if it carries a large fraction, typically 95%, of the jet energy and said to be non-isolated 
otherwise. A comparison between this measured rate and the calculated rate up to O(o) 
Q using the same clustering approach to define the photon yielded a first determination 
of the quark-to-photon fragmentation function accurate at this order. Furthermore, the 
insertion of this measured function into the 0(a) 'isolated' rate, defined as the 'photon' 
+1 jet rate for z > 0.95, yielded a good agreement with the ALEPH data. 

In this fixed order framework, the distribution of electromagnetic energy within the 
photon jet of photon + 1 jet events, for a single quark of charge e q , at 0{a) can be written 

-^ = 2^ 7 (z,M + ^n ( 7 0) Wlogfe) -*)+..., (1) 

a Q dz 7r q ' \ A*f / 

where . . . represents terms which are well behaved as z — ► 1. In the Durham jet algorithm 
and at large z, s cut ~ sz(l — z) 2 /(l + z) ~ p% |4j] where is the transverse momentum of 
the photon with respect to the cluster. The non-perturbative fragmentation function is an 
exact solution at 0{ot) of the evolution equation in the factorization scale 

D q ^(z^ F ) = ^P^(z) log + ZW^/io). (2) 

In this equation, all unknown long-distance effects are related to the behaviour of D q ^(z, fj,o), 
the initial value of this fragmentation function which has been fitted to the data at some 
initial scale hq in [§]. As D q ^(z, p,p) is exact, this solution does not take the commonly 
implemented || resummations of log(/i^) into account and when used to evaluate the photon 
+1 jet rate at 0(a) yields a factorization scale independent prediction for the cross section. 

In the conventional approach, a resummation of the logarithms of the factorization scale 
is performed to all orders in a s and the solution of the evolution equation for D q ^ is 
proportional to log (np/nl). For z < 1, we find that ji 2 F ~ s cut and fip fj^ . The 'direct' 
contribution in eq. ([I]) is therefore suppressed relative to the fragmentation contribution. 
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Figure 1: Parton level subprocesses contributing to the photon + 1 jet rate at 0(aa s ). 



The conventional assignment of a power of l/a 8 to the fragmentation function is in this 
case clearly motivated, this contribution is indeed more significant. However, as z — > 1, we 
see that the transverse size of the photon jet cluster decreases such that s cut — > 0. The 
hierarchy s cut ~ fi 2 F and fi 2 F ^> /Xq is no longer preserved and both contributions in eq. ([[]) 
are important. Large logarithms of (1 — z) become the most dominant contributions. Being 
primarily interested in the high z region, in M it was chosen not to impose the conventional 
prejudice and resum the logarithms of fip a prioribut to work within a fixed order framework, 
to isolate the relevant large logarithms. A detailed comparison of the evaluation of the photon 
production cross section in the conventional and fixed order formalisms can be found in ||. 

We have performed the calculation of the 0(a s ) corrections to the 'photon' + 1 jet rate 
using the same democratic procedure to define the photon as in [@, [§]. The details of this 
fixed order calculation have been presented in [0. We shall here limit ourselves to outline 
the main characteristics of this calculation, to summarize the results and to show how these 
compare with the available experimental data. 



2 The calculation of the 'photon' + 1 jet rate at 0(aa s ) 

The 'photon' + 1 jet rate in e + e~ annihilation at 0(aa s ) receives contributions from five 
parton-level subprocesses displayed in Fig. [|. Although the 'photon' + 1 jet cross section is 
finite at 0(aa s ), all these contributions contain divergences (when the photon and/or the 
gluon are collinear with one of the quarks, when the gluon is soft or since the bare quark-to- 
photon fragmentation function contains infinite counter terms). All these divergences have 
to be isolated and cancelled analytically before the 'photon' + 1 jet cross section can be 



2 



evaluated numerically. 



Within each singular region which we have defined using a theoretical criterion s m i n , 
the matrix elements are approximated and the unresolved variables analytically integrated. 
The evaluation of the singular contributions associated with the process 7* — > qqg'j is of 
particular interest as it contains various ingredients which could directly be applied to the 
calculation of jet observables at next-to-next-to-leading order. Indeed, besides the contribu- 
tions arising when one final state gluon is collinear or soft, there are also contributions where 
two of the final state partons are theoretically unresolved. The three different double unre- 
solved contributions which occur in this calculation are: the triple collinear contributions, 
arising when the photon and the gluon are simultaneously collinear to one of the quarks, the 
soft/ collinear contributions arising when the photon is collinear to one of the quarks while 
the gluon is soft and the double single collinear contributions, resulting when the photon is 
collinear to one of the quarks while the gluon is collinear to the other. A detailed derivation 
of each of these singular real contributions and of the singular contributions arising in the 
processes depicted in Fig. |I](b)-(d) has been presented in M. 

Combining all unresolved contributions present in the processes shown in Fig. |l|(a)-(d) 
yields a result that still contains single and double poles in e. These pole terms are however 
proportional to the universal next-to-leading order splitting function Pffl || and a convolu- 
tion of two lowest order splitting functions, (-P g y ®Pj^ )■ Hence, they can be factorized into 
the next-to-leading order counterterm of the bare quark-to photon fragmentation function 
present in the contribution depicted in Fig. |l](e), yielding a finite and factorization scale 
{Hf) dependent result fTj. 

We have then chosen to evaluate the remaining finite contributions numerically using the 
hybrid subtraction method, a generalization of the phase space slicing procedure [ID, |TT[]. This 



latter procedure turns out to be inappropriate when more than one particle is potentially 
unresolved. Indeed, in our calculation we found areas in the four parton phase space which 
belong simultaneously to two different single collinear regions. Those areas cannot be treated 
correctly within the phase space slicing procedure. Within the hybrid subtraction method 
developed in flPSf , a parton resolution criterion s min is used to separate the phase space into 
different resolved and unresolved regions, but, rather than assuming that the approximated 
matrix elements are exact in the singular regions, the difference between the full matrix 
element and its approximation is evaluated numerically in all unresolved regions. The non- 
singular contributions are calculated using Monte Carlo methods like within the phase space 
slicing scheme. 

The numerical program finally evaluating the 'photon' + 1 jet rate at 0(aa s ) contains 
four separate contributions. Each of them depends logarithmically (in fact as \og 3 (y m - m )) on 
the theoretical resolution parameter y min = s m \ n /Q 2 . The physical 'photon' + 1 jet cross 
section, which is the sum of all four contributions, must of course be independent of the 
choice of ymm- In Fig. ^ we see that the cross section approaches (within numerical errors) 
a constant value provided that y m { n is chosen small enough, indicating that a complete 
cancellation of all powers of log(y min ) takes place. This provides a strong check on the 
correctness of our results and on the consistency of our approach. 
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Figure 2: 0(aa s ) individual contributions (left) and sum of all 0(aa s ) contributions (right) 
to the photon + 1 jet rate for a single quark of charge e q such that ae 2 q = 2tt, a s (N 2 — l)/2N = 
2n using the Durham jet algorithm with y cut = 0.1, and integrated for z > 0.7. 



3 Results 

A comparison between the measured 'photon' + 1 jet rate [@ and our calculation yielded 
a first determination of the quark-to-photon fragmentation function accurate up to 0(aa s ) 



13| . This function, which parameterizes the perturbatively incalculable long-distance effects, 
has to satisfy a perturbative evolution equation in the factorization scale fip. Indeed, the 
next-to-leading order fragmentation function can be expressed as an exact solution of the 
evolution equation up to 0(aa s ) 0, 

2l /,,2\ 
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[z) <8> D(z,fi ) +D(z,no). (3) 



The initial function D(z, yU, ) has been fitted to the ALEPH 1 jet data @ for for the 



jet resolution parameter y cnt = 0.06 and a s (M 2 ) = 0.124 to yield [], 

2 

Z^^/io) = ^ (-P g ( 7 0) (^) log(l - zf + 20.8 (1 - z) - 11.07) , (4) 
where fi = 0.64 GeV. The next-to-leading order (MS) quark-to-photon fragmentation func- 



tion (for a quark of unit charge) at a factorization scale /ip = Mz were shown in |13| and 



1 Note that the logarithmic term proportional to P$ (z) is introduced to ensure that the predicted z 
distribution is well behaved as z — > 1 101 . 
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Figure 3: The integrated photon + 1 jet rate above z = 0.95 as function of y cut , compared 
with the 0(a) and 0(aa s ) order calculations including the appropriate quark-to-photon 
fragmentation functions. 



compared with the lowest order fragmentation function obtained in f2|. A large difference 
between the leading and next-to-leading order quark-to-photon fragmentation functions was 
observed only for z close to 1, indicating the presence of large log(l — z). 

Moreover, a comparison between the ALEPH data and the results of the 0(aa s ) calcula- 
tion using the fitted next-to-leading order fragmentation function for different values of y cut 



can be found in |7|, |L| . The next-to-leading order corrections were found to be moderate for 
all values of y CVLt demonstrating the perturbative stability of our fixed order approach. To 
test the generality of our results, we have considered two further applications: the 'isolated' 
photon rate and the inclusive photon distribution which we shall now briefly present. 

Using the results of the calculation of the photon + 1 jet rate at (D(aa s ) and the fitted 
quark-to-photon fragmentation function, we have determined the isolated rate defined as the 
photon + 1 jet rate for z > 0.95 in the democratic approach. The result of this calculation 
compared with data from ALEPH and the leading order calculation JIJ is shown in 
Fig. |3|. It can clearly be seen that inclusion of the next-to-leading order corrections improves 
the agreement between data and theory over the whole range of y cnt . It is also apparent 
that the next-to-leading order corrections to the isolated photon + 1 jet rate obtained in 
this democratic clustering approach are of reasonable size indicating a good perturbative 
stability of this isolated photon definition. 
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Figure 4: The inclusive photon energy distribution normalized to the hadronic cross section 
as measured by the OPAL Collaboration compared with the 0(a) and (D(aa s ) order calcu- 
lations including the appropriate quark-to-photon fragmentation functions determined using 
the ALEPH photon + 1 jet data. 



The OPAL collaboration has recently measured the inclusive photon distribution for final 
state photons with energies between 10 and 42 GeV [14||. They have compared their results 



with the model estimates of |j[ [15| and found reasonable agreement for /ip ~ Mz in all cases. 
Fig. |] shows our (scale independent) predictions for the inclusive photon energy distribution 
at both leading and next-to-leading order. We see good agreement with the data, even 
though the phase space relevant for the OPAL data far exceeds that used to determine the 
fragmentation functions from the ALEPH photon + 1 jet data M. 



4 Conclusions 



In summary, we have outlined the main features of the calculation [0 of the 'photon' + 1 jet 
rate at 0(aa s ). Although only next-to-leading order in perturbation theory, this calculation 
contains several ingredients appropriate to the calculation of jet observables at next-to-next- 
to-leading order. In particular, it requires to generalize the phase space slicing method 
of [|K], [□]] to take into account contributions where more than one theoretically unresolved 
particle is present in the final state. The 'photon' + 1 jet rate has then been evaluated for a 
democratic clustering algorithm with a Monte Carlo program using the hybrid subtraction 
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method of fll2| . The results of our calculation, when compared to the data on the 
'photon' + 1 jet rate obtained by ALEPH, enabled a first determination of the process 
independent quark-to-photon fragmentation function at 0(aa s ) in a fixed order approach. 
As a first application, we have used this function to calculate the 'isolated' photon + 1 jet 
rate in a democratic clustering approach at next-to-leading order. The inclusion of the QCD 
corrections improves the agreement between theoretical prediction and experimental data. 
Moreover, it was shown that these corrections are moderate, demonstrating the perturbative 
stability of this particular isolated photon definition. Finally, we have computed the inclusive 
photon energy distribution and found good agreement with the recent OPAL data. 
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